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statement  of  Problem  Studied 


The  problem  studied  is  the  search  for  a common  logical  approach 
to  elastic  and  plastic  behavior.  In  particular  the  work  involved 
development  of  the  formxalation  of  a set  of  rate  type  constitutive 
equations  which  exhibits  both  elastic  response  and  plastic  yield 
as  well  as  transition  between  the  two  types  of  behavior,  both  in 
loading  and  unloading.  The  relation  to  classical  linear  elastic- 
perfectly  plastic  behavior  was  studied  and  strain  hardening  was 
considered.  Formulation  of  the  hydraulic  autofrettage  problem 
in  generalized  plane  strain  was  studied.  Plane  strain,  as  well  as 
open  and  closed  end  conditions  were  considered.  Computation  for 
the  open  end  condition  was  begun  for  the  incompressible  case. 
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Sxjmmary  of  Scientific  Work 

Introductory  Note:  We  summarize  here  briefly  the  work  in 

References  [1]  and  [2],  (Reports  2 and  5).  In  addition, 
we  report  in  some  detail  the  progress  on  the  problem  of 
hydraulic  autofrettage  during  the  duration  of  the  project. 

1.  Introduction.  The  purpose  of  the  efforts  on  this  project  was 
to  formulate  properly  a unified  set  of  constitutive  equations 
for  elastic  and  plastic  behavior  and  to  study  methods  of  applying 
the  formulation  to  advantage  in  solving  problems  of  autofrettage. 
During  the  duration  of  the  project,  a formulation  was  brought 
to  completion  and  published  [l].  Methods  were  advanced  for  solving 
the  hydraulic  autofrettage  problem  in  generalized  plane  strain 
as  a step  toward  developing  methods  of  applying  the  equations. 
Computations  were  carried  out  to  test  the  method. 

Althoxigh  the  incompressible  plane  strain  is  considered  too 
simple  to  be  of  great  interest  in  itself,  some  study  was  made  of 
it  initially  in  order  to  test  some  aspects  of  writing  computer 
programs  for  the  equations.  Next  the  cases  of  open  end  and  closed 
end  conditions  were  taken  up  for  generalized  plane  strain. 
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Using  to  advantage  the  lanified  nature  of  the  equations,  we  were 
able  through  mathematical  analysis  to  reduce  these  problems  to 
initial  value  problems  for  a set  of  ordinary  differential  equations, 
which  then  can  be  easily  and  cheaply  solved  by  an  automatic 
digital  computer.  A computer  program  was  written  and  computations 
were  carried  out.  The  computations  demonstrated  clearly  that 
the  method  does  give  an  automatic  transition  to  yield  during 
loading  and  an  automatic  return  to  elastic  behavior  during  un- 
loading. This  is  where  the  work  stood  at  the  time  the  project 
ended.  Although  we  report  numerical  results  which  show  the 
transition  to  and  from  yield,  the  program  written  needs  further 
refinement  before  the  numbers  obtained  are  presentable  as  repre- 
senting the  actual  predicted  stress  distribution.  These  refine- 
ments are  straightforward,  involving  improved  integration 
techniques  and  better  error  bounds,  and  will  be  carried  out  in 
the  future.  The  short  duration  of  the  project  did  not  allow  time 
for  these  refinements.  The  results  are,  however,  successful  in 
what  they  were  intended  to  test,  as  will  be  presented  below. 

2.  Development  of  the  Unified  Equations.  The  results  of  the 
investigations  on  the  unified  equations  was  published  in  the 
International  Journal  of  Engineering  Science  under  the  title  of 
’•A  Unified  Theory  of  Elasticity  and  Plasticity*'  [l].  A pre- 
sentation was  also  made  at  the  fourteenth  annual  meeting  of  the 
Society  of  Engineering  Science  xjnder  the  title  of  "Unified 
Constitutive  Equations  for  Elastic  and  Plastic  Behavior"  and 
an  extended  abstract  was  published  in  the  proceedings  of  that 
meeting  [2].  We  shall  present  here  a brief  summary. 
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A rate-theory  equation  was  obtained  which  relates  stress- 


r 


deviator  s. . and  strain  deviator  C. ..  The  equation  has  the  form 

Pi  1 r 1 p 


s 


id 


^(P) 


;^(P)  ;^'(p) 


where 


and  ^ is  an  unspecified  function.  It  was  shown  by  construction 
that  it  is  possible  to  choose  the  function  so  as  to  make 
equation  (l)  give  rise  to  both  elastic  and  plastic  response  as 
well  as  transition  to  plastic  yield  during  loading  and  from 
plastic  yield  during  unloading  with  a von  Mises  yield  condition. 
Indeed,  yield  corresponds  to  failure  of  a Lipshitz  condition 
with  its  possibility  of  multiple  solutions  and  elastic  behavior 
results  from  viniqueness  of  solutions  when  the  Lipshitz  condition 
holds.  During  loading,  which  is  characterized  by  positive  shear 
work  rate,  yield  is  highly  stable.  During  unloading  (negative 
shear  work  rate)  yield  is  unstable.  Thus  the  assumption  that 
the  stable  solution  is  followed  implies  a return  to  an  elastic 
regime  upon  the  onset  of  unloading  from  yield.  (This  notion  is 
readily  incorporated  into  computer  solutions. ) During  yield, 
the  basic  constitutive  equations  automatically  become  the  Prandtl- 
Reuss  equations.  Furthermore  it  was  shown  how  ^ could  be  chosen 
so  that  the  elastic  response  approximates  to  any  desired  degree 
a linear  elastic  behavior.  An  extension  to  strain  hardening 
was  constructed. 
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3*  The  Problem  of  Hydraulic  Autofrettage.  It  must  be  landerstood 
that  the  work  at  this  stage  is  intended  to  develop  new  methods 
of  solving  autofrettage  problems.  It  is  natural,  then,  first  to 
attack  old  problems  by  a new  method  before  launching  on  new  ones 
in  order  to  test  the  method.  Such  has  been  the  intent  of  the  work 
done  so  far.  We  have  during  the  duration  of  the  project  brought 
the  work  to  the  point  of  establishing  the  feasibility  and  ease  of 
applying  computer  techniques  to  the  minified  equations.  In 
particular,  there  are  three  matters  which  have  been  tested: 

1)  advantages  of  the  unified  equations  in  analyzing  a problem 
in  preparation  for  computation,  2)  the  transition  to  yield  dviring 
loading  in  a computer  solution  and  3)  the  use  of  the  stability 
criterion  to  determine  when  there  is  persistence  of  yield  and 
when  there  is  a return  to  elastic  behavior.  We  feel  that  these 
tests  have  been  successfully  met. 

Work  was  done  on  the  case  of  generalized  plane  strain  with 
different  end  conditions.  Analysis  for  both  the  incompressible 
and  compressible  cases  were  performed.  However,  to  date  time 
has  permitted  only  the  incompressible  case  to  be  brought  to  the 
stage  of  computation.  We  shall  now  discuss  generalized  plane 
strain. 

The  (physical)  components  of  displacement  u^^  in  cylindrical 
coordinates  r,  6,  z are  taken  to  be 

U = Uj,,  V = Uq,  w = u^, 

and  for  generalized  plane  strain  we  assume 

u = u(r),  V = 0,  w = Az + B. 

For  the  incompressible  case,  the  equation  of  incompressi- 
bility gives 
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(2) 


whence 


„ C Ar 
^ = r-  T 


(3) 


where  C is  a constant  of  Integration.  The  stress  cr.  . and  the 

J 

stress  deviator  s. ^ are  related  by 

O'.  . = -'D  O . . + S . . 


s . . =0. 
11 


\diere 

P = - 

Solutions  are  then  sought  for  which  the  off-diagonal  components 
of  Sj^^  vanish,  and  we  write  more  simply  s^  for  Sq  for  Sqq 


and  s^  for  s^^' 


Equilibrium  then  gives 
-dp  9sp  s - Sq 

+ — i + -£ 2 = 0 

dr  dr  r 


(4) 


or 


Bo;,  s - 
^ ^ = 0. 
or  r 


(5) 


We  took,  as  in  reference  [l],  5^  to  satisfy 

_ 1 

- 4(i-  ^ 


where  n is  an  integer,  ^ is  a constant  giving  initial  shear 
modulus  and  k is  a constant  giving  shear  yield  stress.  The 
quantities  were  normalized  or  non-dimensionalized  in  such  a way 

p 

as  to  lead  to  the  same  result  as  taking  2p  = 1 and  2k  = 1. 

(This  can  be  done  by  replacing  s^^  by  s^^/k  and  u^  by  |iu^/k. ) 

We  shall  then  assume  that  the  quantities  actually  represent  their 
non-dimensionalized  equivalents,  whence  we  write 


6 


1 

^ ' = (1  - ^ 

and  obtain  ^ as  in  reference  [l].  Then  equations  (l)  and  (3) 
give  for  the  incompressible  plane  strain 


s 


z 


^ A*  1 


(6) 


Let  a be  the  inner  radius  and  let  b be  the  outer  radius  of 
the  cylinder  which  is  to  be  stressed.  For  all  cases  we  shall 
have  oj,  = 0 at  r = b.  Let  P be  the  hydraulic  pressure  of  the 
fluid  inside  the  cylinder.  Then  for  all  cases  we  require  that 
o^  = - P at  r = a.  The  conditions  distinguishing  the  three  cases 
are  respectively  the  following: 

Case  I:  Plane  strain:  A = 0. 

Case  II:  Open  end  condition: 

1 r dr  = 0 (?) 

'^a 

Case  III:  Closed  end  condition: 
f ^ 2 

\ r o-^,  dr  = (8) 

J a 

We  now  discuss  the  cases  more  specifically: 

Case  I:  Plane  strain,  A = 0.  In  this  case  we  can  take 
= 0 and  the  problem  reduces  to  solving  the  simple  equation 
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r 

with  = -P/i/5  for  which  we  could  take  C as  a parameter: 

dp  ^ I 


Equation  (9)  was  solved  numerically  at  a number  of  r locations 
using  a Hamming  Predictor-Corrector  Scheme  [3]-  The  following 
instruction  was  added:  If  P i 1,  replace  P by  N,  where  N was 
taken  to  be  a number  slightly  less  than  unity,  and  n was  given 
several  integer  values,  1,  20,  32,  64.  N was  given  values  .99 
and  .9999.  Achievement  of  yield  occurred  as  expected  during 
loading,  and  a return  to  elastic  behavior  occurred  during  un- 
loading. Indeed,  the  only  difficulty  which  was  experienced 
occurred  for  the  physically  unimportant  small  n = 1,  where  it 
was  found  that  for  the  values  of  N chosen,  extra  accuracy  would 


have  to  be  introduced  into  the  integration  scheme  in  order  to  get 
a transition  to  the  proper  elastic  regime.  This  phenomenon  was 
ascribed  to  the  relatively  slow  change  in  slope  near  yield  for 
n = 1 and  did  not  occur  for  the  higher  values  of  n.  This 
observation  leads  to  some  interesting  mathematical  questions, 
which  were  not  pursued  at  this  stage  since  they  did  not  affect  the 
cases  of  interest  for  which  the  scheme  worked  well.  Once  P is 
known  as  a function  of  r and  C,  we  have  s^  = -p//5,  s^  = 0,  Sq  = 
P//?,  and  P can  be  determined  by  solving  (8)  with  o^(b)  = 0 and 
using  P = -cy^(a).  The  incompressible  plane  strain  case,  then, 
is  basically  one  dimensional  and  quite  trivial.  We  go  on  to  the 
other  two  cases. 
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Cases  II  and  III:  Open  and  Closed  End  Conditions;  Note  that 


O'  = O'  + s - s 
2 r 2 r 


whence  we  can  write  for  the  open  end  condition  (7) 
b b 

j ro;,  dr  + ( r(s^-  s^)dr  = 0. 


'a 


For  the  closed  end  condition  (8),  we  get  similarly 
.b  .b 


rcr;  dr  + | r(s2-s^)dr  = 


a"P 


(10) 


(11) 


Integration  by  parts  and  the  conditions  at  the  inner  and  the 
outer  radius  give 

(12) 


^ Hv.  Pa 
r o^  dr  = 


'a 


Substitution  of  dcr^dr  from  (5)  into  (12),  replacement  of 
Sq  by  -s^-  s,  and  substitution  of  the  result  into  (10)  and  into 
(11)  gives 


I 


rs,dr.-£|- 


a 


for  the  open  end  condition  and 


(13) 


1 r s^  dr  = 0 
a 

for  the  closed  end  condition  respectively. 

Integration  of  (5)  and  use  of  the  inner  and  outer  radius 
conditions  gives 
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f*’ 


dr. 


(14) 


For  the  open  end  condition,  we  obtain  by  differentiating  (13) 
and  (14)  and  substituting  (6) 


2* 

-a^P 


= TA  + LC 


-P=  MA  + QC 


where 


p 

' ■ !.'  [? 


?■ 


dr 


(15) 


M = 


3 s s + ^ s ^ 
^ r 2 2 2 


dr 


Q = - 


D 

1 r-2p , 1 (,/r~^  p ’ll 

7Li^vr'^  -y)\ 


2 2 
4 s^  + 4 s„s^  + s„ 
r r 2 2 


J dr. 


(16) 


For  the  closed  end  condition  we  obtain  similarly 
0 = TA+  LC 
-P  = MA  + QC  . 

For  the  open  end  condition  we  can  solve  (15)  for  A and  C in 
terms  of  P,  substitute  the  result  into  (6)  and  write,  after 
dividing  both  sides  by  P 
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1 P T 1 


^ j^^(  2s2  + s^s^ ) + ^ s^s^ 


1 fpL  1 [/<  ;:2n  pIFt/^. 

"ar'KCTEj-p^^ p 


r®z 


(18) 


The  problem  of  generalized  plane  strain  for  either  case, 
open  or  closed  ends,  can  now  be  turned  into  an  initial  value 
problem  with  ordinary  differential  equations;  We  discretize 
the  interval  [a,b]  by  subdividing  it  into  a number  of  equal 
intervals.  In  practice  we  used  100  Intervals.  Let  the  points 


of  subdivision  be  r^,  1 = 0,...,100  with  a = rQ,  b = Let 

the  values  of  s^,  and  s^  at  the  mesh  points  r^^  be  s^  and  s^ 


11 


rr^ 


respectively,  i = 0,...,100.  We  then  express  all  integrals 
appearing  in  (16)  in  terms  of  the  mesh  values  s^,  s^  through 
a numerical  quadrature  formula.  (In  practice,  the  trapezoidal 
rule  was  used. ) This  allows,  then,  either  the  equations  (17) 
(open  end)  or  the  equations  (18)  (closed  end)  to  be  replaced  by 
a set  of  202  ordinary  differential  equations  for  which  we  take 
as  initial  conditions  s^  = s^  = 0 at  P = 0.  The  bovindarv  and 
end  conditions  have  been  absorbed  and  the  problem  is  reduced  to 
an  initial  value  problem.  Such  problems  are  easily  auid  cheaply 
solved  by  digital  computers. 

Recall  now  that  in  the  discussion  above  on  the  plane  strain 
case  (case  I),  which  reduced  to  the  one  dimensional  problem  (9), 
there  was  a discussion  about  how  to  treat  P when  it  reached  unity 
in  order  to  avoid  an  error  signal  when  the  numerical  procedure 
approximated  p by  a value  greater  than  one  and  in  order  to  assure 
the  slight  perturbation  required  to  insure  a return  to  an  elastic 
regime  at  the  onset  of  unloading.  For  the  two  cases  now  under 
discussion,  namely  the  open  and  closed  end  conditions,  this 
instmiction  has  to  be  generalized  to  one  which  will  hold  when 
the  problem  no  longer  reduces  to  a non  dimensional  one.  The 
instruction  devised,  then,  was  the  following:  If  p ^ 1,  replace 
®ij  Ns^j/P,  where  N is  a number  slightly  less  than  unity. 
During  loading,  this  perturbation  should  be  dying  out. 

Calciilations  were  carried  out  on  the  open  end  case  in  order 
to  test  the  procedure.  The  results  of  these  calculations  left 
no  doubt  that  the  method  works  well  and  cheaply.  (A  typical  run 
on  the  computer  cost  between  five  and  ten  dollars. ) Yield 
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occurred  during  loading  and  a return  to  elastic  behavior  occurred 

at  the  onset  of  unloading.  Results  are  shown  in  figures  1 and  2, 

in  which  the  ratio  of  the  inner  radius  to  the  outer  radius  was 

taken  to  be  two.  Figure  1 shows  P versus  r during  loading  for 

several  values  of  hydraulic  pressure  P (non  diraensionalized) , 

which  was  brought  up  to  a maximum  value  of  0.7^  in  steps  of 

0.01.  One  can  see  the  yield  region  (p  = 1)  progressing  from  the 

inner  radius  outward  for  the  curves  P = .70  and  P = .74.  Figure 

2 

2 shows  the  behavior  of  P during  unloading:  Note  a return  locally 
to  elastic  regimes  as  P drops  below  its  maximum.  Comparison 
of  figures  1 and  2 shows  that  unloading  behavior  is  quite  different 
from  loading  behavior.  This  is  because  material  points  which  have 
experienced  yield  return  to  a different  elastic  regime  than  they 
were  in  before  yield. 

The  program  needs  refinement  before  the  actual  numerical 
values  of  stress,  in  particular  residual  stress  can  be  reported. 

A simple  Euler  method  was  used,  N was  taken  to  be  .99  and  the 
function  Y crudely  calculated.  There  is  no  difficulty  at 
all  in  making  the  necessary  refinements,  and  they  will  be  done 
in  order  to  control  error  buildup.  The  project  ended  before  they 
could  be  accomplished  and  so  we  report  the  state  of  affairs  at 
that  time. 

We  conclude  now  that  the  essential  nature  of  the  procedure 
has  been  tested  and  shown  to  work.  It  is  easy  to  handle  and 
economical  computation-wise  because  of  the  analysis  that  the 
unified  equations  allowed.  We  feel  that  it  is  extremely 
promising. 
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We  shall  end  with  a brief  discussion  of  the  compressible 
case,  for  which  analysis  has  been  done.  The  analysis  is  similar 
to  that  for  the  incompressible  case  and  we  shall  merely  report 
how  it  differs  from  that  already  presented. 

For  the  compressible  case  we  donnot  have  the  equation  of 
incompressibility  (4).  Instead  we  have 


P = (19) 
where  K is  a constant,  namely  the  bulk  modulus.  Using  (19)  and 
(5),  we  obtain  from  the  equation  of  equilibrium  (4) 


u = - 


Ar  r 

~T  ' W 


i: 


3^,(4  ) - Sq(<J) 


di, 


^ 1 4s2(<^)d4  - Ik  • 

a 


(20) 


The  expression  (20)  for  u is  then  substituted  into  the 

constitutive  equations  (l),  where  again  the  off  diagonal  components 

of  s. . are  assumed  to  vanish.  The  result  is  a set  of  two  volterra 
-*•  J 

integral  equations  of  the  second  kind  for  s^,  and  s^.  Such 
equations  are  classical  and  portend  no  difficulty  in  solving  for 
Sj,  and  s^.  The  resulting  equations  then  replace  (6)  and  thence- 
forth the  procedure  is  the  same  as  in  the  Incompressible  case. 

The  general  philosophy  behind  these  procedures  is  that  the 
larger  the  part  of  the  task  relegated  to  analysis  and  the  smaller 
the  part  to  computation, the  better,  more  accurate  and  cheaper  will 
be  the  computation.  In  other  words,  the  extra  thought  that  goes 
into  analysis  pays  off.  We  are  happy  that  we  have  been  supported 
in  bringing  our  method  to  the  stage  at  which  it  now  stands.  And 
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from  the  evidence  accumulated  thus  far,  we  are  confident  that 
it  will  prove  itself  capable  of  giving  excellent  results  as  it  is 
tested  on  problems  and  applied  to  new  problems.  It  should 
become  an  excellent  tool  for  problems  of  autofrettage  as  well 
as  plasticity  in  general. 
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Figure  Captions 
2 

Loading  behavior  of  p = s. .s. . versus  radial  coordinate 

1 J 1 J 

(r-a)/(b-a)  for  inner  radius  for  b/a  = 2 at  various 
values  of  internal  gas  pressure  P.  This  is  the  incom- 
pressible case  of  generalized  plane  strain  with  open 
end  conditions  normalized  on  yield  stress  so  that 
2k^  =1,  n = 20. 

2 

Same  as  figure  1,  except  that  it  shows  behavior  of  P 
during  imloading. 


Figure  1. 


Figure  2. 


15 


